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ABSTRACT 


An automatic, electronic controller was designed for use with 
adiabatic calorimeters. The controller employed reset action, and 
its behavior at different system gains was analyzed using the 
analytic techniques of feedback control theory. In accord with the 
analysis satisfactory performance was obtained with a calorimeter 
characterized by large thermal lags, but marginal performance was 
obtained with another calorimeter characterized by a rapid thermal 
response, A second controller, which employs reset and proportional 
action gave satisfactory results during preliminary testing for both 


calorimeters. 
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I. INTRODUCTION 


When thermodynamic measurements are made in an adiabatic calori- 
meter, the heat exchange between the sample and its enclosure should 
be zero. To accomplish this goal, it is necessary to observe and con- 
trol the temperature difference between the Sample and its enclosure, 
or shield. 

A calorimetric experiment begins with the sample at a known 
temperature and in thermal equilibrium with its surroundings, the 
shield, A measured amount of energy is then added to the sample. The 
shield control system (controller) must act to re-establish thermal 
equilibrium between the sample and the shield within a reasonable 
length of time and in such a way that little, or no, heat is lost from 
the sample. After equilibrium is established, the temperature of the 
sample is again recorded. During equilibirum periods, the controller 
must supply enough heat to the shield to offset losses by conduction 
through its supporting structure and by radiation to the calorimeter 
jacket, which is maintained at some constant lower temperature by its 
contact with a thermal reServoir or constant temperature bath. During 
each period when the sample is heated, the shield requires additional 
heat in order to increase its temperature at a rate to match that of 
the sample. 

In a typical system, a differential thermocouple is utilized to 
sense the sample~shield temperature error, thus providing information 
by which automatic regulation of the shield temperature is accomplished. 
The error signal is thus inherently a low level signal, which implies 


that the control system have high gain. To maximize sensitivity, it 


is of importance to avoid electromechanical components and the inherent 
deadband which they induce into the control system and to maximize the 
gain, consistent with stability, with which the error signal is intro- 
duced into the control a cemae The most convenient method of supply- 
ing and controlling the shield heat is electrical, such as by electroni- 
cally controlling a current to a resistance heater on the shield. 
Automatic thermal control devices fall into two categories: 
(1) intermittently acting devices, such as a thermostat, and (2) con- 
tinuously acting, or modulation, devices incorporated into a closed 
loop feedback system. The modulating type control action is preferred 
over the on-off action since the feedback path continually senses the 
effect™of the manipulated variable, the shield current, on the equili- 
brium process, Electronic control action is commonly categorized into 
three processes, proportional, rate and reset, any combination of 
which may be incorporated into the design of a modulated shield heater 
current, control sie Proportional control produces a heater cur- 
rent directly related to the magnitude of the input signal from the 
thermocouple. Automatic reset control is an integrating action which 
provides a floating correction proportional to the time duration as 
well as the magnitude of the error from the equilibrium set point. Rate 
control applies a correction proportional to the rate of change of the 


error signal and represents a differentiating action. 


panic. Elwood R. and Moiles, E. D., "Electronic Systems for Temper- 


ature Measurement and Control", Temperature - Its Measurement and Con- 


Pood mine oc iencewandmindusitny evs ll lop. / /2-//35. C.sHerzteld, edeaNew 
York, Reinhold Publishing Corp., (1962). 


“Thid., p. 771-772. 


This paper describes the design and operation of an all electronic, 
dc shield temperature regulator operating in the reset mode. The regu- 
lator operates as one element of a closed loop feedback system whose 
dynamic behavior is investigated using the analytical techniques of 
feedback control theory. The regulator is to be used in an adiabatic 
calorimeter employed in taking specific heat and electrocaloric measure- 
ments on ferroelectric crystals in the neighborhood of the ferro- 


electric transition. 


IL. CONTROLLER DESIGN AND PERFORMANCE 


Since there is interaction between the various parts of the system, 
such as the shield, controller and the sample, the system must be 
treated as a whole in order to design a satisfactory controller. A 
simple, yet effective, tool for visualizing the behavior of a system 
is the block diagram. The blocks represent components, or component 
assemblies, whose physical behavior may be described mathematically by 
algebraic functions under static conditions and by integrodifferential 
equations under dynamic conditions, Figure 1 displays the pertinent 
variables and the interconnecting information streams of a closed loop 
temperature control system. Any difference between the set point 
Signal and the output of the thermocouple is presented to the control 
system as an error Signal. The output of the control system regulates 
the shield heater so as to correct for disturbances. The effect of 
this control action is sampled by the thermocouple sensor and fed-back 
for comparison with the set point, thus completing the feedback loop. 

Since an equilibrium state for the system requires no control action 
other than the maintenance of a steady state value of current to the 
shield heater, which corresponds to a null or Zero signal from the 
thermocouple, the two-input system may be reduced to a simpler one- 
input system, and rearranged as shown in Fig. 2. For proper control 
action, the error signal must be presented to the controller with the 
proper sign, as indicated in the diagram, so that the loop is repre- 


sentative of a negative feedback system. For example, if the 
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FIGURE 1. General Closed Loop Control System 
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FIGURE 2. Single Input Closed Loop Control System with 
Laplace Notation 
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disturbance was such that the sample has just béeen eee ee the polarity 
of the thermocouple signal must activate the controller to increase the 
heating current supplied to the shield heater so that the shield temper- 
ature may increase to match that of the sample. 

For the controller discussed in this thesis, reset action was 
chosen as the method of control, instead of utilizing a combination of 
two or more types of control action, such as proportional plus rate or 
the more conventional use of proportional action by itself. An inte- 
grating type controller was chosen since a practical integrating device 
could be constructed simply and economically using solid state devices, 
and because the thermal lags associated with the typical adiabatic 
shield allows stability in the feedback loop to be achieved by a lag 
control network, The control action is accomplished by storing charge 
on a capacitor and modifying the charge strored by means of the error 
signal. The voltage appearing at the integrator output drives a cur- 
rent amplifier with the resistance heater on the shield for a load, 

Several factors relative to the interaction between the controller 
and the calorimeter shield must be considered in order to obtain 
proper reset action. While a large integrator time constant permits 
distrubances of a long duration to be integrated without saturation of 
the integrator, a large time constant may inhibit control action, re- 
sulting in overdamping. A decrease of the time constant will result 
in increased hunting, and may also lead to loss of control. Preampli- 


fication of the error signal increases the sensitivity of the controller, 


Breas of the sample and the measurement of its temperature are 
accomplished by separate, independent circuity which is external to 
the thermal equilibrium feedback Loop. 


UZ 


but too large a value of preamplification may cause increased hunting 
if the heater current is changed too quickly and too much. If thermal 
lags associated with the diffusion of heat through the shield are 
large enough, additional damping may not be necessary in the controller; 
while if they are not, damping may be added to the control system in 
the form of a compensation network. If system damping is larger than 
desired, the compensation network will take the form of a lead net- 
work to lessen the damping. If thermal lags are too large, the con- 
troller cannot be made to give a rapid response, 

Thus, a simple reset action controller, composed of a preamplifier, 
a lead or lag compensation network, an integrator, and a current amplLi- 
fier, behaves in a quantitatively complicated fashion. In designing 
the controller, the techniques of automatic control theory, to which 
the block diagram has provided a brief introduction, provide some in- 
Sight into a qualitative design solution for the performance of the 
temperature controller in a system whose thermal constants are known, 
Final values of controller components may be selected after suitable 
operating experience has been obtained using the engineering estimates 
arrived at as a result of the control theory analysis, 

It is easier to manipulate the Laplace transforms of the blocks in 
Fig. 2 than it is to work in the time domain. Thus, the following 
analysis uses the rules of transform algebra found in any standard 


ae , : 
EouENOneautomatic control theory. 7 The two basic assumptions in the 


See R. N., Introduction to Automatic Control Systems, Wiley, 
boo 


ae Valkenburg, M. E., Network Analysis, Prentice-Hall, 1955. 
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analysis are (1) the system elements can be described by lumped-para- 
meters, as opposed to distributed parameters, and hence the system may 
be described by ordinary Linear differential equations rather than 
partial differential equations whose transforms are more complex, and 
(2) the response functions of system elements are Linear. Non-linear 
functions are Linearized by assuming the function can be approximated 
by just the first-order term of its Taylor series expansion, With 
reference to Fig. 2, one can define the following quantities in trans- 
ron noOtatLon: 


G(s) is the thermocouple transfer function; 
H(s) is the controller transfer function; 
E(s) is the Laplace transform of the error, ¢€(t); 


B(s) is the Laplace transform of the primary feedback signal, 
bie); 
C(s) is the Laplace transform of the actuating signal, c(t)and 


Q(s) is the Laplace transform of the input disturbance, q(t). 


These variables are related by 


E(s) = G(s) x C(s), 
B(s) = H(s) x ECs), 
C(s) = Q(s) - BCs). 


These equations may be combined to eliminate B(s) and C(s), 


Eis) _ G(s) = ya) 
Q(s) 1 - G(s)H(s) 


where W(s) is the system transfer function. 
The thermocouple has a first order time constant response, since it 
is basically an energy storage element in series with a resistance to 


flow. The forward path transform has the following form: 
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The controller tramstorm, HaseGefig J @2e41s Construct camusanmenc 
transforms of the individual controller stages - preamplifier, compen- 
sation network, integrator and buffer/amplifier. A high gain preampli- 
fier with adjustable gain, K,, has a small time delay, L/e, ants 


response so that its transform is 


ee 
A ste 


A first-order phase lead compensation network has the form 





Wisi aioe 2 


An ideal operational amplifier, which is used to construct an analog 
integrator, has gain, input impedance, and bandwidth which approach 
infinity and zero output impedance: Although real operational 
amplifiers have Limitations, which complicate the form of the inte- 
grator transform, currently available, solid state, de operational 
amplifiers have characteristics which approach the ideal so that the 


complete integrator transform may be approximated by a simple function, 


GL sw 7 L/RC 


The buffer and current amplifier stages following the integrator have 
transforms which are pure gains t2king the forms K, and it respect- 
ively. If a stage does not load the previous stage, the complete 


controller transform can be written 


K(s_ + aje 


7 Cee ean ~ s(s + b) 


K.K.K /RG sec! 
Bp 


2a 
I 


~ 
I! 


with A ; 


aalepey.. G., "Analog Integration", Instruments end Control Systems, 
v. 42, p 49-50, January 1969. 
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The system transfer function may now be specified as 


—_ Can c (stb) (ste)s 


~ 4 
L+ GH st 4 (btdte)S” + (betbdtde)S” + (bde + Kce)S + Keea 


where, the open loop gain is given by 


> a Kec (sta) 
oor S(st+b) (s+d) (ste) : 


and 1 + GH is defined to be the system characteristic equation. Figure 
3 is a root locus plot of the characteristic equation which shows the 
loci of the roots in the s plane as the gain, K, increases from Zero 
toward infinity. For K = 0 all the roots lie in the stable, left-half 
of the s plane. As K is increased toward infinity, the poles move to- 
ward the zeros (there are 3 zeros at minus infinity) along the dipole 
branches and the pole pair near the origin enters the complex plane as 
a conjugate pair. The imaginary axis is a region of marginal stability 
and the system will exhibit "ringing" at a critical value for K which 
places poles on this axis. From Fig. 3, one sees that oscillationgyig? 
occur sooner for large K in a system with small thermal constants, c 
and d, than in a system in which the shield-thermocouple arrangement 
has a large thermal lag. 

The system behavior is examined by applying a disturbance, Q, in 


the form of a step signal, T/s, and examining the resultant response, 


DeiCsiecwili(Seae) 
— 4 3 2 
S + (btd+e)S + (bcthdtde)S + (bde + kce)s + Kcea 


The general procedure is to factor the denominator into four roots, 
expand E using partial fractions, solve for the residues and perform 


the inverse transform to the time domain with the aid of a table of 
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Laplace transforms. The solution of the quartic with particular numer- 
ical coefficients is straightforward, albeit Pet ous. but it is too 
difficult to attempt for the general case. However, certain approxi- 
mations may be made which reveal the nature of' the solutions for two 
cases; K small, say K~ 10, and large K in the region where ringing 
occurs. 

For small K (and small s), the pole due to the preamplifier lag is 
considered ‘remote’ by virtue of its location in the remote part of the 
s plane. The residue at this pole is small compared to the other 
residues and, in addition, the time constant of this pole is very short 
when compared to the other time constans so that its contribution to 
the system behavior is negligibly small. For small K, with G, ~ K 


A A 


and Q = T/s, the error behaves like 


nvr 


ey I (a4 3 al aye mf exp[-t(b + d)/3}. 


At the output of the termocouple, this equation has an initial value of 


q (0) = K e(0) = T(2b - d)/3. 


The disturbance is converged to zero, and since q(o) =T, by defini- 


tion, there exists a requirement that 


Beet det eon /2. 
For large values of K, at large s, the pole of the pole-zero dipole 
branch is nearly cancelled by the zero and the residue of this pole 
becomes negligibly small. In this case, the compensation network has 


no effect, and the system behavior is characterized primarily by the 


l prslagieon, R. S., Handbook of Mathematical Tables and Formulas, 
Boe Joe Handbeokebublishers, 19572 
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remote pole and the conjugate pole pair. Shunting of the integrator 
input by a larger b/a ratio prevents extension of the compensation net- 
work to very high gains. No attempt was made to derive the general 
solution for e(t); however, the error does converge to zero since the 
behavior of the system is the sum of a dying exponential from the re- 
mote pole and a sine term multiplied by a dying exponential from the 
conjugate pair. This behavior is valid until K is large enough to 
move the pole pair across the imaginary axis into the region of in- 
Seaprlity. An application of the Routh Criterion for stabil ome 

shows the critical gain to be 

d 
Cc 


K x E 


As K is increased beyond the critical value, the system nonlinearities 
assume a dominate role. At a high enough value of K, the system re- 
Sponse will ultimately be bounded by a limit cycle operation describ- 
able by rapid shutoff. Thus, automatic shield control is accomplished 
through a combination of continuous and intermittent control as K is 
increased through several orders of magnitude. 

Using the results of the previous analysis, construction of the 
controller procedes in a straight-forward manner. For measuring and 
amplifying the output of the differential thermocouple, which is 
typically fractions of a microvolt, 2 Keithley model 149 millimicro- 
voltmeter was chosen. The model 149 is highly stable, has nanovolt 
sensitivity and contains built-in zero suppression to cancel Spurious 


potentials generated along the thermocouple transmission line and at 


the amplifier input. The amplifier output drives the controller shown 


en Valkenburg, pervect wpe SL? 


Le, 


schematically in Fig. 4. As the discharge of the capacitor will 
normally be controlled by the input impedance of the buffer, it is 
then required that the leakage resistance of the capacitor exceed the 
input impedance of the buffer. This requirement can be met by a 
variety of capacitors. Since a large capacitance is necessary to pro- 
duce a reasonable time constant, a tantalum capacitor was used. The 
capacitor may not be reduced with reduction compensated by an increase 
aot Ry» as the input current offset, drift and noise of the operational 

ee ”, 
amplifier Al are proportional to R, and these effects may become 
troublesome enough to require separate compensation if R, is further 
NO Nile 
increased, The output stage can deliver 4W of continuous power to 
a load drawing 140 ma, when a heat sink is used with Q.° DS) is a 10S% 
temperature sensitive device which stabilizes Q, during sustained 

eel 

operations at higher currents. The controller was mounted on a 
small etched circuit plug board. A standard 5 in x 19 in rack panel 
contains the controller, the necessary power supplies, an output cur- 
rent monitor, and a battery source for manually driving the controller. 


The input is Zener Limited to a + 10.2 volts to protect Al from 


over Loading. 


ine A 741, Fairchild Semiconductor, 313 Fairchild Drive, 
Mountain View, Cadi fommecay, 


Bede. C. V., "Review of Operational Amplifier Principles," 
Application Bulletin - July 1966, Fairchild Instrumentation, 475 


Ellis Street, Mountain View, California 


nssaet en manufactured by Texas Instruments Inc., Post Office 
Box 5012, Dallas, Texas 
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The controller constants may be computed from the component values 
in Fig. 4. The output stage voltage gain, Kp, is derived using an 


approximate two parameter equivalent current. 


tes EST 


SS hl + +h oR, ~ by /e, FR, 7 18 for Ry = 208 
ie te- 4e 


The transistor parameters were measured with a Tektronix Model 57/5 


Curve Tracer, and found to be he. =~ 17.5K and hh. w 1/5. The buffer 


fe 
stage voltage gain expression is derived from the JFET equivalent 
; ols: , : ; 
Cir Cue and the gain Ka» is computed using the transistor parameters 
given by the manufacturer's data sheet, 
RS 
SOF 


+ Rite + 1) 


= ' — 
; 0.9 where R! R_//{R, i eC aes h, RA]. 


The Keithley time constant was measured by observing the output volt- 
age rise time to a step input on the 100 UV scale with a Tektronix 


model 454 CRO and found to be 
1/é@ ~eeensec 


For the compensation network used 


s -1 
a = 1/R,C. = 1 sec : and b= 1/(R,//R,) = 1,2 sec 


If the thermal constants for a particular system are known, the 
idealized behavior for an assumed input may be predicted from E = QW, 


An examination of the orders of magnitude of the various terms in the 


earealimane J. and Halkias, C. C., Electronic Devices and Circuits, 
pp. 340-341, McGraw-Hill, 1967. 


meatal pp. 404-405. 
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transfer function reveals which element(s) of the controller dominates 
its behavior for a given gain. The critical gain for stability of the 
system can only be calculated if the thermal delays of the calori- 


meter are specified. 





FIGURE 5. Controller Etched Circuit Board 
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The controller was tested with two calorimeters of very different 
thermal design. The first calorimeter, for which the most extensive 
testing was carried out, was for use near 120°K and was characterized 
by sluggish thermal response. The second, a developmental calori- 
meter, was designed for use near room temperature and was character- 
ized by very small thermal time constants. 

The lower temperature calorimeter contained a cylindrical Cu 
shield, with dimensions 2.0 in o.d. x 0.016 in wall x Fails in length, 
and mass of 98.7 gms. The shield was isolated from the calorimeter 
jacket by a cyclindrical Cu furnace with dimensions, 2h78 in’ id ies 
0.002 in wall x i in length, and mass of 38 ems. Both cylinders 
screwed onto a 47 em Gu lid which was silver brazed to a % in i.d. 0.01 
in stainless-steel tube 8.5 cm long. The entire assembly was sealed 
within a vacuum jacket which, in turn, was immersed in a bath of Liquid 
nitrogen. The heater consisted on an 80 2 manganin coil around the 
lid stem in series with a 27 (Q) manganin coil around the outside sur- 
face of the furnace cylinder so that most of the heat appeared at the 
lid. Current distribution between the two heater coils was externally 
controlled by a potentiometer to trim gradients and the entire heater 
assembly appeared as a 200 Q load at the controller terminals. The 
thin walled stainless-steel tube was characterized by a thermal 


resistance, 


24 


: C ee meer: ; 
where Kss = 0.17 W/cm - K. At Tf 120 K radiation was relatively 


14,15 


( 


all 
The thermal lag of the calorimeter appears in the system transfer 
function of the previcus chapter through G, the thermocouple transfer 
function, which is d=teérmin=3 cxpzrimentally. The testing technique 
measures the transient response of the thermocouple by comparing the 
measured temperature reeding to 2a theoretical ramp change in tempera- 
ture produced by imposing 4 step chenge in Current through the shield 
fe Le 16 
heater under -diabatic conditions. Neglecting the IR drop at the 
thermocouple-shield junction, the calorimeter tested had time con- 


stants given 2pproximarely bv 
L/w 6 sec and tI/d w~3 sec. 


These constants zre for a 0.002 in diam chrcemél-alumel thermocouple 
bonded to the shield with GE7031 varnish. When these values are in- 
serted into the previous th-acry for e(t), the maximum allowable gain 


far aontinuews zontrol action is 


Since this thermcvouple produc-s an average of 15uV per degree of 





Et eae W. and May, L. & YS. ea. G2. 5 SZ (1967). 


cae W. end May, L. F., Leberstory Notebock, January 1966 - 
June 1967, Naval Pestersadust- School, Monterey, California 93940, 
unpub lished, 

16 . i ee , 
Green, S. J. and Hunt, i. W., "“Sccurb®y and Response ef the 
Thermocouples tor Surface and Fluit Tomperature Measur=nents," 
Temperature-Its Me=sursment and Control in Scienc= anc industry, 
pp. /08-/11. 


2) 


differential, the Keithley should be operated on the 1 UV Scale where 
Ky = tee for maximum sensitivity. At this value of Ky» full scale 
ringing occurred and the controller entered a limit cycle. However, 
the performance of the controller on the 1 pV scale, as measured by 
the thermal response of the sample, was acceptable as shown in Fig. 6 
and was an improvement over a proportional controller which was pre- 
viously used with this see 

For too large a Sample heating step, corresponding to coarsely 
spaced data points, the controller applied maximum available power 
initiating an extended Limit cycle observed as an overshoot on the 
sample recording. The maximum sample heating step was experimentally 
determined and depended on several factors; for example, the current 
gain of the output stage, the heat capacity of the sample and the 
proximity of the thermocouple junction to the sample heater. 

Drift of the sample temperature for the sample used to obtain Fig. 
6 at a point about 5°K below the ferroelectric anomaly was improved 
from 0.045°K/hr to 0.006°K/hr, measured over a 9 hour period, after 
a few minutes of careful adjustment of the Keithley zero suppress. 

The controller was less satisfactory for the second calorimetric 
system. The second system was characterized by a light mass Al con- 
struction and coupling to the bath was by a short brass Link with Rint 
8.3 K/W. In addition, radiation coupling of the shield to the bath 
was much more important (wT?) as was the increased coupling of the 


sample to the shield. The constants, c and d, were not determined for 


this calorimeter but the time response was at least a factor of 10 
Ly 
Reese and May, Phys. Rev. 162, 612 (1967). 
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*The vertical scale represents temperature increase and the horizonta] 
scale shows clapsed time with 50 sec between divisions 


FIGURE. 6. Sample Trace of Temperature Versus Time Illustrating 
Equilibrium of Potassium Dihydrogen Phosphate for a 
.025°K Step Beginning at T = 118.15 1°K. 
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faster. A sample recording exhibited unsatisfactory oscillations 
tracable to the ringing of the thermocouple signal. The trace oscilla- 
tions were minimized, but not entirely eliminated, by reducing the inte- 
erator time constnat by a factor of 10, Limiting the maximum current 
delivered to the load and reducing the gain, Kye Despite these changes 


the system still performed only marginally; however, the system could 


LO overcome this problem, a proportional control stage was added 
and a second controller was assembled which operated in both the reset 
and proportional modes. The test circuit is shown in Fig. 7 with 
approximate component values. This network has exhibited extremely 
good stability with the Cu calorimeter. Control action was so posi- 
tive that the thermocouple error never exceeded 2 UV during sustained 
Sample heating under maximum power. The performance of this controller 
in the second system was an improvement over the reset action con- 
troller; however, the sensitivity of the calorimeter still exceeded 
the ability of the controller to eliminate hunting on the sample trace 
at the required gain level of Ky = 107 under all conditions, although 
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stable operation was obtained with Ky Seo xe Cae 
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